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Abstract 

Density functional calculations are performed for ground [He]2s^ ^S"^, and three metastable bound 
excited states, Is2s2p2 ^pe^ Is2p3 ^S", Is2s2p3p ^P^ of Li" and [He]2s2p2 ^pe^ [He]2p3 ^S", Is2s2p3 
of Be~ each. The work-function-based exchange potential is used, while the correlation effects 
are included by employing the Lee- Yang-Parr potential. The relevant nonrelativistic KS equation 
is solved by means of a generalized pseudospectral discretization scheme offering nonuniform and 
optimal spatial grid. Computed total energies, radial densities, selected density moments, as well 
as two transition wavelengths (Is2s2p2 5pe^is2p3 ^S" of Li", [He]2s2p2 ^p^^ [He]2p3 "^S" of Be") 
show reasonably good agreement with the available theoretical and experimental data. The term 
energies show an absolute deviation of 0.007-0.171% with the largest deviation being observed for 
the even-parity ^P state of Li^. The transition wavelengths of Li^, Be~ are calculated within 
0.891 and 0.438% of the experimental values. This offers a simple practical route towards accurate 
reliable calculation of excited states of anions within density functional theory. 
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I. INTRODUCTION 



Atomic negative ions are representative of weakly-bound (typically having binding ener- 
gies one order of magnitude smaller than the neutral atoms or positive ions) fragile quantum 
systems. Due to the weak central field experienced by its outer valence electron, correla- 
tion is dominant and can no longer be treated as a small perturbation to the independent 
particle picture. Thus faithful estimation of binding energies, fine structure or excitation 
energies of negative ions pose considerable challenges to theoreticians and experimentalists. 
For example, very good precision on total energy is necessary to obtain modest values of the 



as well as experimental works have 
2I for review) during the past three 



electron affinity. An impressive amount of theoretica_ 
been reported in the literature (see for example, 
decades and interest in this area continues to grow further. 

Our focus is on the ground and excited states of negative ions of two most extensively 
studied systems (after possibly H~ and He~) from alkali and alkaline-earth groups, viz., Li, 
Be. Ever since the three long-lived bound states of Li~ were reported through configuration- 
interaction (CI) calculations 3|, significant attention has been paid on their characteriza- 
tion and understanding. These three are respectively: the ground state [He]2s^ ^S*^ and two 
core-excited high-spin metastable bound states, Is2s2p^ ^P^, ls2p^ ^S°, each lying below the 
corresponding Is2s2p and ls2p^ ^P parent states of Li. Thereafter a host of calculations 
including multiconfiguration Hartree-Fock (MCHF) [s^, multiconfiguration interaction {g] as 
well as the recent saddle-point variational Monte Carlo (VMC) [8], and other methods 
have been employed to determine the upper bound, relativistic and nonrelativistic energy, 
transition wavelength, oscillator strength, lifetime, fine structure, hyperfine parameter, etc., 
of these states. In general, these findings show good agreement with the experimental obser- 
vations jsl, [lo|. In contrast to Li^ ground state, addition of a 2p electron to the ground state 
of Be ls^2s^ yields an unbound Be" [He]2s^2p ^P° ground state; however, three metastable 
states [He]2s2p^ ^P^, [He]2p^ and Is2s2p^ were predicted, of which the first two have 



been observed experimentally UHlSj. Theoretical calculations include a variety of meth 



ods, e.g., superposition of configurations 1^, CI 



15 



16|, state-specific theory [l7|, MCHF 
19|, full core plus correlation [20|, Rayleigh-Ritz variation method [21|], VMC [8], etc., 
among others. 

In the last four decades, density functional theory (DFT) [22| has emerged as one of the 
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most powerful and promising tools for electronic structure and dynamical studies of many- 
electron atoms, molecules, solids. While it enjoyed remarkable success for the ground states, 
same for excited states came much later. In the last decade, a work-function-based DFT 



prescription has 
for example, 23 



Deen shown to be fairly successful for general atomic excited states (see 



24| and the references therein). This proposed the use of a nonvariational 



local work- function-based exchange potential [25j], which is computationally advantageous 



compared to the nonlocal HF potential. Usua 



relation functional of Lee- Yang-Parr (LYP) 26| was employed in conjunction. Recently this 



ly, the gradient and Laplacian-included cor- 



was extended for higher excitations through a generalized pseudospectral (GPS) implemen- 
tation which allowed to solve the resultant radial KS equation in a non-uniform, optimal 
spatial grid accurately efficiently j23|-j3li. This enabled us to study, for the first time, 
within the framework of DFT, a broad range of important physical processes in atomic ex- 
cited states such as multiple excitations, valence as well as core excitations, satellite states, 
ligh-lying Rydberg states etc. However with the exception of 8 triply excited states of He~ 
23[, all of these works dealt with either neutral atom and cations. The current work aims 
at investigating the spectra of two most important and widely studied case of negative ions 
through the above DFT method, with an objective to assess and judge the performance of 
this simple single-determinantal approach in the context of excited states of negative ions. 
In other words can a density-based approach as ours, produce satisfactory and physically 
meaningful results for these intricate excited states; and if so how good do they perform in 
comparison to the traditional wave function-based methodologies, which are prevalent for 
their studies? To answer these questions, we investigate the term energies, radial densities, 
transition wavelengths as well as expectation values of some position operators for 4 states 
of Li", VIZ., [He]2s2 Is2s2p2 sp^^ is2p3 5S°, Is2s2p3p 5p^ and 3 states of Be", viz., 
[He]2s2p^ ^P*^, [He]2p'^ and Is2s2p'^ ^S°. Comparison with available theoretical and ex- 
perimental results are made, wherever possible. The article is organized as follows: Section 
II gives basic elements of the formalism as well as its numerical implementation. Section III 
makes a discussion on the results while we end with a few concluding remarks in Section IV. 
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II. METHOD OF CALCULATION 



The methodology as well as the numerical solution of the desired radial Kohn-Sham (KS) 
equation through GPS discretization scheme has been presented before [23,24,27-31]; hence 
not repeated. Only essentials details are given. Unless otherwise mentioned, atomic units 
employed throughout the article. 

Our interest is in the following single-particle nonrelativistic KS equation, 

1. 



ipiir) = Eiifjiir) (1) 



where three terms in the left-hand side signify respectively the kinetic, electrostatic and 
exchange-correlation (XC) energy contributions, fes(r) contains the nuclear-attraction and 
classical internuclear Coulomb repulsion as, 

-e.(r) = -- + / ^dr' (2) 
r J I r — r' I 

where Z corresponds to the nuclear charge. The local nonvariational exchange potential 
relates to the work required to move an electron against the electric field Sx{r) arising out 



of its own Fermi-hole charge distribution, Px{r, r'), and given by the line integral 25 ] 



V:Ar 



rS,{r)-dl (3) 

J oo 



e,ir) = / ^-<V;'"''„; dr. (4) 



where 

J \Y — T\ 

For well-defined potentials, work done must be path-independent (irrotational), which is 
rigorously satisfied for spherically symmetric systems. This potential can be calculated ac- 
curately as the Fermi hole is known exactly in terms of the single-particle orbitals. Working 

within the central-field approximation, ^j(r) = Rniir) Yimi^t), employing a suitable corre- 

I i 

lation functional (here we use LYP potential, [26;]), solution of the KS equation produces a 
self-consistent set of orbitals, which gives the electron density as, 

p(r) = Ei^«wr- 

i 

The radial KS equation is solved accurately and efficiently by means of the GPS method, 
which has shown remarkable success for structure and dynamics of Coulombic singular sys- 
tems such as atoms, molecules as well as other difficult and stronger singularities such as 



TABLE I: Calculated total energies of ground and excited states of Li^, Be^ along with the 
literature data. Numbers in the parentheses denote absolute per cent deviations. See text for 
details. 



Ion 


State 






-E(a.u.) 




X-only 






XC 


This work 


Ref. 


This work 


Ref. 


Li- 


[Hc]2s2 iS-^ 


7.4278(0.005) 


7.4282" 


7.4984(0.009) 


7.4553\7.5008'', 7.4991"' 




Is2s2p2 -^P^ 


5.3640(0.006) 


5.3643"= 


5.3925(0.171) 


5.3866''' = , 5.3833^5.3865^ 5.3863'' 




ls2p^ 5go 


5.2223(0.004) 


5.2225"= 


5.2608(0.137) 


5.2561''' = , 5. 2536^5. 2560A 5.2558* 




Is2s2p3p 


5.3289 




5.3683(0.007) 


5.3679= 


Be- 


[He]2s2p2 4p= 


14.5078(0.008) 


14.50909 


14.5806(0.062) 


14.5779'', 14.5716^ 14.57089, 14.5769^ 




[He]2p3 4s° 


14.3272(0.002) 


14.32759 


14.4081(0.049) 


14.4063'', 14.4010^ 14.40029 




Is2s2p=* •^S" 


10.4279(0.009) 


10.42889 


10.4758(0.092) 


10.4662<=, 10.46159, 10.4711' 



'Ref. [si]. ''Ref. [l^. '^Ref. Q]. <^Ref. "Ref. Q. '^Ref. il. SRef. ''Ref. 'Ref- ^Ref- Q- ''Ref. fl. 



the spiked harmonic oscillators, Hulthen and Yukawa potentials etc. 27|]- 3l|. In most of 
these cases, this has either outperformed the best results published so far for those systems 
or offered results of comparable accuracy to the best calculations. This is achieved through 
the following steps: (i) approximate a function by an Nth order polynomial (we use Leg- 
endre) fN{x) in such a way that the approximation is exact at collocation points xj, i.e., 
fN{xj) = f{xj) (ii) a transformation r = r{x) employed to map the semi- infinite domain 
r G [0, oo] onto the finite domain xG [— 1, 1] (iii) introducing an algebraic nonlinear mapping 
r = r{x) = L{l + x)/(l — x + a), with a = 2L/rjnax, followed by a (iv) symmetrization proce- 
dure. Finally this leads to a symmetric eigenvalue problem which can be easily diagonalized 
efficiently by using standard available libraries (like NAG) to generate accurate eigenvalues 
and eigenf unctions, for both low as well as higher levels. Another notable feature is that 
this nonlinear optimal discretization maintains similar accuracy at both smaller and larger 
distances with significantly smaller number of points (all calculations done with 250 radial 
points), than those needed in usual finite-difference or finite element methods but at the 
same time promises faster convergence. A convergence criterion of 10~^ and 10^^ a.u. was 
imposed for the potential and energy to obtain all the results presented here. By performing 
a series of test calculations, a consistent set of GPS parameters were chosen (a=25, N=250, 
'^max = 200 a.u.) which produced "stable" converged results. 



5 



III. RESULTS AND DISCUSSION 



Table I reports our calculated density functional term energies for ground and excited 
states of Li~ and Be~. We consider 4 states for the former, viz., [He]2s^ ^S^, Is2s2p^ ^P*^, 
[He]2p3 53°, Is2s2p3p ^P^; and 3 for the latter, [He]2s2p2 ^pe^ [He]2p3 ^S", Is2s2p3 ^S". 
For each of these states, two sets of energies are given; viz., X-only (noncorrelated) and 
XC (correlated). All of these states have been quite extensively studied, except the core- 
excited even-parity Is2s2p3p ^P of Li~, which has been reported only lately. The latter 
is also the only case for which HF result is unavailable. Existing literature data are given 
for comparison. To put results in proper perspective, in parentheses the respective per 
cent deviations are given; for X-only case these are relative to the lone literature results 
in column 4; for XC case, these are with respect to the recent variational Monte Carlo 
(VMC) [^] values except for the 4th states of Li^, where such a result does not exist; and 
this is given in reference to the recent saddle-point calculation of In VMC approach 
the authors employed explicitly correlated wave functions having a Jastraw factor and a 
multideterminant model wave function to account for the dynamical and nondynamical 
correlation effects respectively. Present X-only result of ground state of Li~ is higher from 
the accurate HF calculation of 32(] by only 0.0004 a.u. Our XC energy value is in fairly good 



agreement (slightly above) with the accurate correlated results reported through MCHF-n 

1 I 

expansion considering all expansions [32], as well as the VMC method [8]. However the 
earlier result of [l^ seems to be in considerable disagreement with these two and ours. The 
present X-only results for the next core-excited high-spin even-parity ^P"^ and odd-parity 
states of Li~ again show excellent agreement with the recent HF estimates of 811, while the 

n n 

XC energies match well with the VMC J8|], extensive CI [3i], variational multiconfiguration 
calculation jo], saddle-point j^, MCHF [5], etc. Note that the present XC energies for these 
two states are lower than all of these reference results by 0.171 and 0.137%, respectively 
and constitute the two instances giving maximum deviations in our calculation. Since our 
X-only results are virtually of HF quality, this overestimation is presumably caused by the 
correlation potential employed. It is worth noting here that even though a KS equation is 
solved with the work-function exchange and LYP correlation potential, the procedure is not 



subject to a variational bound 



23 



241. We note that some of these correlated calculations 

ji 

are highly elaborate and extensive; for example, [6[ used a 45 angular component 1004-term 
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wave function, j3| used a 320-term CI, [t] used 7-50 angular spin components and 541-1298 
linear parameters for the former, and it is gratifying that current results exhibit a rather 
small and acceptable discrepancy. The next even-parity ^P^ state has been investigated only 
recently [tI, and our energy shows very good matching with theirs. They used a saddle-point 
restricted variational method with accurate multiconfigurational wave functions built from 
STO basis sets. As already mentioned, Be~ has no [He]2s^2p ground state; and three 
metastable bound states found in the discrete spectrum are given in this table. [He]2s2p^ 
^PS [He]2p3 and Is2s2p3 63° states lie below the Be Is22s2p ^P, Is22p2 3p and ls2sV 
^P excited states. The X-only energies again matching excellently with the HF results 15 |. 
A large number of accurate and sophisticated theoretical results exist for the correlated case. 
For example CI calculation of jl5] included single, double, triple and quadrupole subshell 
excitations. The first two states have also been studied through a method of full core plus 
correlation and restricted variation approach 20| . Our XC energies show better agreements 
with literature results in this case than those for Li~; however once again falling below the 
reference values for all the three states as in Li^. Also a combined Rayleigh-Ritz and method 
of restricted variation result exists for the state 21|. 

Next, two calculated transition wavelengths (in nm) of Li^ Is2s2p^ ^P'^ — t- ls2p^ and 
Be^ [He]2s2p^ ^P^ — )■ [He]2p^ '^S° are reported in Table II. A variety of theoretical and ex- 
perimental results have been published, which are given for comparison. Note that the latter 
includes relativistic effects, whereas present calculation is nonrelativistic. Generally there 
seems to be good agreement for these quantities among the literature values. The present 
results are lower than the experimental values of 9] and [l^ by 0.891% and 0.438% respec- 
tively for Li~ and Be~. In parentheses we give the corresponding wavelengths calculated 
from the X-only results. As expected they show considerable difference. 

Now Fig. 1(a) depicts the single-particle radial density plots for Li~ ground state [He]2s^ 
^S"^, and two excited states Is2s2p^ ^P^, ls2p^ ^S°. Close to the nucleus, the latter two 
states making a family of similar radial densities, possess considerably smaller charge density 
than the ground state, mainly due to the presence of two core Is electrons in the ground 
state and one Is electron in the latter two. However this scenario reverses in the region 
enclosed between the first minimum and up to r = 4 a.u., where the ^P^ and states 
show larger values than the ^S"^ state. The two excited states also appear to separate out 
in this intermediate region with ^P*^ showing larger peak value than the ^S°. After that, at 
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TABLE II: Comparison 


of transition 


wavelengths for Li , Be with 


literature data. 1 


a u =27 2113834 eV hc= 


197 3269602 eV 


nm. X-only values are enclosed, in 


Daren theses 






Wavelength (nm) 




Transition 


This work 


Other theory 


Experiment 


Li- Is2s2p2 5pe _5. is2p3 5s° 


345.96 


346.06'',349.12*,349.0^348.98'' 


349.07^349.0^ 




(321.55) 






Be- [He] 2s2p2 4p<= _j> [He] 2p3 


■*S° 264.14 


267.19,265.4'',265.370',265.32J, 265.04* 


265.30l',265.318™,265.331" 




(252.29) 








r (a.u.) r(a.u.) 

FIG. 1: Radial density plots for several states of (a) Li^ and (b) Be~. 

large distances, all these three states show similar behavior with the ground state having 
greater charge density and decaying rather slowly than the two excited states; ^P^ dying 
out fastest. Similar radial density plots for all three excited states of Be~ are given in Fig. 
1(b). In the vicinity of nucleus, the ^P^, states show similar behavior and have much 
larger charge densities compared to the state, again presumably because the former two 
have two core Is electrons while the latter has only one. At the intermediate distance after 
the first minimum and up to r = 3.5 a.u., this situation changes with having the largest 
charge density. Also the ^S°, ^P^ states branch out in this region with the latter showing 
slightly larger peak value than the former. After that at larger distances, all the three states 
decay in a similar pattern with ^S" having higher values and oozing out slowly than the 
other two and decaying out first. These behaviors in electron density are also reflected 
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TABLE III: Calculated expectation values in a.u. Numbers in the parentheses are taken from 



Ion State (r-2) (r"!) (r> (r^) 

Li- [He]2s2 iS'= 30.458(30.19) 5.920(5.8977) 11.846(11.718) 73.713(71.30) 

Is2s2p2 Spe 19.022(19.04) 4.173(4.1548) 11.933(11.7410) 53.197(57.726) 

Is2p3 5S° 18.325(18.23) 4.094(4.0475) 12.148(12.9592) 70.686(72.73) 

Be- [He]2s2p2 56.927(56.7) 8.622(8.5525) 11.131(12.478) 62.927(66.86) 

[He]2p3 4s° 56.083(55.7) 8.485(8.4140) 12.173(14.135) 82.928(89.97) 

Is2s2p3 34.468(34.6) 6.276(6.2137) 10.038(10.5831) 30.512(34.785) 



in the selective radial expectation moments (r^^), (r^) and (r^) values of these 6 

states of Li~, Be" as compiled in Table III. Some of these quantities have been recently 
reported in the VMC calculation of js], quoted in parentheses for comparison. Generally 
there is good agreement between our results with them, especially for the negative moments 
and discrepancy tends to increase for the positive moments. However the general trend is 
similar as theirs. The negative expectation values of Li~ are reduced for both the excited 
states (^P*^, ^S*^) compared to the ground state, and having closer values to each other, due 
to the close resemblance of their corresponding radial densities close to the nucleus. ^P'^ 
state with lowest energy showing larger (r~^) and (r^^) than the ^S°. The positive moments 
however correspond more to the charge densities at intermediate-to-large distances, and show 
complicated behavior. For Be~, the negative moments are smaller for term compared to 
the quartet S and P terms, with the latter two having similar values, as they have the same 
[He] core. As in Li~, the P state (with lower energy) gives higher values of the negative 
moments compared to the S state. 



IV. CONCLUSION 



A density-based formalism was employed for the first time to calculate the state energies, 
radial density as well as selected density moments of three excited states of Li~ and Be~ 
in addition to the ground state of the former. A work-function-based exchange potential 
plus the nonlinear LYP correlation potential was used, while the resulting KS equation was 
numerically solved by an optimal, nonuniform, spatial GPS discretization scheme accurately. 
Additionally two transition wavelengths of Li~, Be~ were computed. Comparison with the 
existing theoretical as well as experimental data were made. For all practical purposes 
our X-only results are of HF quality and with correlation included, they produce excellent 
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agreement with the hterature data. The observed absolute per cent deviation in the term 
energies between our XC calculations and the correlated results in the literature remains 
within 0.007-0.171% only, and is mostly attributable to our current inadequate knowledge of 
the delicate correlation effects in these strongly-correlated systems, besides the fact that the 
presently employed LYP potential was designed primarily for the ground states. Construc- 
tion of the exact many-body correlation functional has been one of the most forefront areas 
of DFT ever since its rejuvenation in the works of Kohn-Sham in 1964 and unfortunately still 
remains one of the most important unresolved issues in DFT. However, even though its exact 
form remains elusive, many accurate forms have been proposed over the decades suitable 
for different processes in atoms, molecules, solids. Note that this is a single-determinantal 
approach and is far more simpler to implement compared to the benchmark sophisticated 
and elaborate CI or multiconfiguration type formalisms within the traditional wave function 
based methodologies; yet as our results indicate, this can lead to quite good-quality results, 
which of course is the appeal and usefulness of DFT-based methods. At this stage, we note 
that, success of such a relatively simpler method on the present and previous atomic excited 



states done before, e.g., in |23|, Ud, l27|, l28|, may encourage one to explore its validity and 
applicability for other physical and chemical systems such as molecules. While this could 
lead to an attractive and direct route towards molecular excited states, as discussed earlier. 



(see for example, 27| and references therein), a straightforward extension is difficult; as it 
relies on the tacit assumption of spherical symmetry in the density, which is not the case 
for molecules. Finally, while the LYP correlation potential employed has provided rather 
excellent results for the present inherently correlated systems and in past studies, neverthe- 
less it is unclear whether the non-dynamical correlation contribution is described by LYP 
properly or for these states such contributions are small. No definitive answer is possible and 
more detailed and sophisticated studies would be required for this. In conclusion, we have 
proposed an accurate and reliable method for the excited states of negative atoms within 
the framework of DFT. 
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